arXiv: 1506.01705v2 [hep-ph] 2Jul2015 


ZU-TH-16/15 
June 2015 


On the breaking of Lepton Flavor Universality in B decays 

Admir Greljo 0 ’ 6 , Gino Isidori a c , David Marzocca a 

(а) Physik-Institut, Universitat. Zurich, CH-8057 Zurich, Switzerland 

(б) Faculty of Science, University of Sarajevo, Zmaja od Bosne 33-35, 

71000 Sarajevo, Bosnia and Herzegovina 
(c) INFN, Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 


Abstract 

In view of recent experimental indications of violations of Lepton Flavor Universal¬ 
ity (LFU) in B decays, we analyze constraints and implications of LFU interactions, 
both using an effective theory approach, and an explicit dynamical model. We show 
that a simple dynamical model based on a SU (2 )l triplet of massive vector bosons, 
coupled predominantly to third generation fermions (both quarks and leptons), can 
significantly improve the description of present data. In particular, the model de¬ 
creases the tension between data and SM predictions concerning: i) the breaking 
of T—fj, universality in B —»• D^Iv decays; ii) the breaking of fi—e universality in 
B —> KI + I~ decays. Indirectly, the model might also decrease the discrepancy be¬ 
tween exclusive and inclusive determinations of \ V c b\ and \ V u b\- The minimal version 
of the model is in tension with ATLAS and CMS direct searches for the new massive 
vectors (decaying into r + r~ pairs), but this tension can be decreased with additional 
non-standard degrees of freedom. Further predictions of the model both at low- and 
high-energies, in view of future high-statistics data, are discussed. 
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1 Introduction 


Recent experimental data in B physics hint toward deviations of Lepton Flavor Universality 
(LFU) in semi-leptonic decays, both in the case of b —> c charged-current transitions, as well as 
in the case of b —> s neutral currents. The statistically most significant results can be summarized 
as follows: 


• 3.8<r + 2.0 a deviation of r/t universality (£ = r, e) in b —> c transitions, encoded by [1-3]: 


R 


T/t 

D* 


R 


r/t 

D 


B(B -4 D*Tu) exp /B(B -4 D*tv ) sm 
B(B -4 DHu) exp /B(B -4 DHv ) SM 
B(B -4 DTv) exp /B(B -4 Dtzz) sm _ 
B(B -4 D£v) exp /B(B -4 Dlv ) SM = 


( 1 ) 

( 2 ) 


• 2.6(7 deviation of [i/e universality in b —> s transitions [4]: 1 


R 


M/e 

K 


B(B -4 K^ + H )exp 
B[B -4 Ke+e-) exp 


g 2 * * * S[l,6]GeV 


0.7451K ± 0.036 . 


(3) 


1 The result in Eqs. (1) and (2) are obtained using B(B -4 D*rv)/B{B -4 D*iv) exp = 0.323 ± 0.021 and 

B(B —> Dtu)/B(B -4 Dlv) ex p = 0.41±0.05 from the average of Babar [1], Belle [2], and LHCb [3], assuming e//r 

universality in b —> civ decays, as indicated by b —> civ data [5] (see Sect. 3.1), together with the theory predictions 

B{B -4 D*tv)/B(B -4 D*fo)sM = 0.252 ± 0.003 [6] and B{B -4 Dtv)/B{B -4 U/m) S m = 0.31 ± 0.02 [7]. The 

SM expectation of R^/ E is \{R ^/ e )sm — 1| < 1% [8] while, by construction, R^J* = R 7 ^ 1 = 1 within the SM. 
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In addition to these LFU ratios, whose deviation from unity would clearly signal physics 
beyond the Standard Model (SM), 13-physics data exhibit other tensions with SM expectations 
in semi-leptonic observables. Most notably, a ~ 3 a deviation from the SM expectation has been 
reported by LHCb [9] in the so-called P 5 differential observable of B —> K* fi + decays [10]. 
Moreover, in charged current transitions there is a long-standing ~ 2.5cr discrepancy in the 
determination of both | V c b\ and \V u b\ from exclusive vs. inclusive semi-leptonic decays [ 11 ]. 

These deviations from the SM have triggered a series of theoretical speculations about possi¬ 
ble New Physics (NP) interpretations, see in particular Ref. [13-28]. Among these recent papers, 
two particularly interesting observations are: i) the proposal of Ref. [18] to explain both R^ c 
and the P 5 anomaly by means of NP coupled dominantly to the third generation of quarks 
and leptons, with a small non-negligible mixing between third and second generations; ii) the 
observation of Ref. [19] that is natural to establish a connection between R^ e and R r Jy if the 
effective four-fermion semi-leptonic operators are build in terms of left-handed doublets. 

Despite this recent progress, a coherent dynamical picture explaining all the anomalies has 
not emerged yet. On the one hand, a significantly improved fit of experimental data can be 
obtained with a specific set of four-fermion operators of the type J q x Jg, where J q and are 
flavor-non-universal left-handed quark and lepton currents [19,28]. On the other hand, even 
within an Effective Field Theory (EFT) approach, it is hard to believe that this set of effective 
operators is the only relevant one in explicit NP models. In particular, explicit NP models should 
face the tight constraints on four-quark and four-lepton operators dictated by meson-antimeson 
mixing, and by the bounds on Lepton Flavor Violation (LFV) and LF non-universality in pure 
leptonic processes. Moreover, the size of the SM modifications in Eqs. (l)-(3) points toward 
relatively light new degrees of freedom, that could well be within the reach (or already excluded) 
by direct searches at the LHC. 

In this paper we present an attempt to build a simplified coherent dynamical model able to 
explain, at least in part, these violations of LFU. The guiding principle of our construction is the 
idea that the J q x J? effective operators are generated by the exchange of one set (or more sets) 
of massive vector bosons that transform as a SU{2)l triplet, and that are coupled to both quark 
and lepton currents. This hypothesis allows us to establish a connection between quark-lepton, 
quark-quark, and lepton-lepton effective operators. We further assume that the flavor structure 
of the new currents is consistent with an approximate U(2) q x U{2)^ flavor symmetry acting on 
the first two generations of quarks and leptons, along the lines of Ref. [29]. 

Under these assumptions we proceed with two main steps: i) we analyze the low-energy 
constraints (and the corresponding phenomenological implications) on the complete set of four- 
fermion operators generated within the model; ii) we discuss the additional constraints due 
to electroweak precisions test and collider searches, following from the specific choice of the 
mediators. 

We find that, after taking into account all the existing constraints, the proposed model can 
still provide a significantly improved fit as far as low-energy observables are concerned. The 
most serious constraint on the model follows from the searches performed by ATLAS and CMS 
on new heavy neutral states ( Z') decaying into t + t~ pairs. However, as we will discuss, the 
tension with direct searches can be decreased with additional non-standard degrees of freedom, 
whose net effect is the enhancement of the Z' decay width and the corresponding suppression 
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of the Z' —> t + t~ branching ratio. The tension can be further reduced in the limit where the 
assumption of narrow resonances (T -C M), that is implicit in all present direct searches, no 
longer holds. 


2 The model 


2.1 Step I: four-fermion operators 

Our main assumption is that all the non-standard four-fermion interactions can be described by 
the following effective Lagrangian 


A£ 


(T) 

4 / 


2 niy 


ja ja 

”[i i 


(4) 


where is a fermion current transforming as a SU(2)l triplet, built in terms of SM quarks 
and lepton fields: 

j; = 9q^ij (Oh^Qi) + geXfj (P Llll T a L{) . (5) 

Here are Hermitian flavor matrices and, by convention, A 33 = A 33 = 1. 

We define Q l L and U L to be the quark and lepton electroweak doublets in the flavor basis 
where down-type quarks and charged-leptons are diagonal. We assume an approximate U{2) q x 
U(2)i flavor symmetry, under which the light generations of Q l L and L 1 l transform as 2 q X 1^ and 
1 q x 2p, respectively, and all other fermions are singlets. We further assume that the underlying 
dynamics responsible for the effective interaction in Eq. (4) involves, in first approximation, 
only third generation SM fermions (the left-handed l q x H fermions). In this limit, the flavor 
couplings in Eq. (5) are \f- = di 3 S 3 j. The corrections to this limit are expected to be generated 
by appropriate U{2) q x U(2)n breaking spurions, connected to the generation of subleading terms 
in the Yukawa couplings for the SM light fermions. 

In the quark case, the leading U{2) q breaking spurion is a doublet, whose flavor structure is 
unambiguously connected to the CKM matrix ( V ) [29]. We can thus expand A?- as follows: 


A l = 5 i3 5 3j + (eife V 3j + e\V^5 3j ) + e 2 (V£%) + . • • , V 3j = V 3j - S 3j V 3j , (6) 

with €2 = 0(e f). As we will discuss below, low-energy flavor-physics data imply e 7 ; -C 1- 

The breaking structure in the lepton sector is less clear, given the intrinsic ambiguity in 
reconstructing the lepton Yukawa couplings under the (natural) assumption that neutrino masses 
are generated by a see-saw mechanism. 2 As we will discuss below, low-energy data are compatible 
with the hypothesis that the leading breaking terms in the lepton sector transform as doublets 
of U(2) t . ‘ 

2 An attempt to build a consistent neutrino mass matrix starting from an approximate U(2)e symmetry broken 
by small U(2)e doublets has been discussed in Ref. [30]. 
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Among the four-fermion operators generated by the model, the ones most relevant to flavor 
phenomenology are: 


ACiV = - 


AC 


(T) 

FCNC 


AC 


(T) 

AF=2 



I a Lll A) +h.c. 
- (y\«v % U\ 


8 rriy i_ 


(VW% (uhA 


a4fv = -ACKMAtWhA) 


AC 


(T) _ 

LFU — 


9e 

8 rriy 


i-^ibKd + ^adKb)(^L'y^L)( i/ L'y^L) ■ 


\i \(oa 


(7) 

L - ) , (8) 

(9) 
( 10 ) 
(ID 


2.2 Step II: simplified dynamical model 

(T) 

In order to generate A C\j in a dynamical way, we introduce the heavy spin-1 triplet, 

(a = 1,2,3), following the general simplified Lagrangian proposed in Ref. [42]. By means of 
this approach we can describe both models in which the new vector is weakly coupled, such as 
gauge extension of the SM, and strongly coupled models, such as Composite Higgs models. The 
simplified Lagrangian reads 

C v = + gnV^H^T a i D, H) + V“j; , (12) 


where T a = a a /2, - D v V* and D tl V°' = + ge abc W b V„ c . * * 3 

By integrating out at the tree-level the heavy spin-1 triplet and keeping only effective oper¬ 
ators of dimension < 6, we obtain the effective Lagrangian 


r*d=6 


2 rriy 


ja ja _ 

^ (L ” fl 


9h 

2 rriy 


lH ] T a i D M H)(H*T a i D, 


H) - ^L(H ] T 
™v 


•fA 

l i D u 


H)J° 


(13) 


(T) 

By construction, the first term is A C\J in Eq. (4). The second term, in the unitary gauge, is 
simply 

- % (i+") 4 ■ < i4 > 

This term induces an unphysical (custodially-invariant) shift in the W- and Z-boson masses, 1 , 
that can be reabsorbed by a redefinition of v, and deviations in the Higgs interactions to W 
and Z bosons. The latter are well within the existing bounds for the relevant set of parameters. 

4 With respect to Ref. [42] we dropped interaction terms with two or more insertions of V£. While such 

terms can be relevant for double production, they do not contribute to the low-energy effective Lagrangian at the 

dimension-6 level and are thus largely unconstrained by low-energy data. 

4 Within the full model of Eq. (12) this corresponds to a mass mixing between the SM EW gauge bosons and 
the heavy vector triplet. The relative shift in the heavy vector masses my is only of 0(gj I m\ V v 2 /my) 
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The last term, instead, describes non-universal deviations in the Z and W couplings to SM 
quarks and leptons that lead to non-trivial constraints on the parameter space of the model (see 
sect. 4.1). 


3 Low-energy implications of the four-fermion operators 


3.1 New physics effects in charged currents 


Since the new interactions are purely left-handed, in the case of charged currents their effect is 
simply an overall (flavor non-universal) rescaling of the SM amplitudes: 


r)b—> c _ A(b —> C iy*)sM+NP 
A(b —> c UF*)sm 


1 + RqX^ 



where 


Ro 


2 M T ) 
9q9l m W = G F 

g 2 rny G|, m 


Using this expression, the LFU breaking ratio in Eq. (1) assumes the form 


(15) 


(16) 


R 


■It¬ 

'D* 


1 + 2R 0 Re 


1 - 



Vcs(\ q bs y + v cd (K d ) 
v cb 


1 + 2R 0 , 


(17) 


where we have assumed lA^gJ <C 1 and |A?-| <C |V^*V 3 j|. The first condition is required by the 
smallness of deviations from the SM in b —> dtv decays (see below), the second condition follows 
by the consistency of the bounds from A F = 2 amplitudes (see Sect. 3.2). We are thus able to 
fix the overall strength of the new effective charged-current interaction (compared to the Fermi 
coupling): 

R 0 = (R$- l) =0.14 ±0.04 . (18) 

The model predicts the same violation of r/t universality for all type of b — > c and b —>• u 
transitions. This implies, in particular, R r J) = R T Jy* ■, that is perfectly consistent with the 
experimental result in Eq. (2). 

In principle, violations of LFU universality are expected also between b —> c{u)gv and b —> 
c(u)eis modes. Unfortunately, it is very difficult to get experimental bounds on the latter using 
published data, since most of the high-statistics semi-leptonic analyses are performed combining 
H and e modes [5]. Using the PDG fit for the combined B±/B° sample [5], that is separated for 
H and e modes, we deduce that deviations between T(b —> c{u)/iv) and T(6 -» c(u)ev) as large 
as ~ 2% are allowed by present data. Within our model, we expect 


T(b —> c{u) jUZ/)sM+NP 
r(6 -»■ c(u) ei/)sM+NP 


l + 2 Ro(A^-A £ ee ) ■ 


(19) 


The strong constraints on LFU involving only quarks and leptons of the first two generations 
(tv and K decays, CKM unitarity, and n decay [31]) implies |Af e | <C |A^ M |. As a result, the 
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( 20 ) 


constraints on /j-e charged-current LFU can be used to set the approximate bound 

«*' 5 °- 07 (^) • 

that justifies having neglected A* ;i in Eq. (17). After this bound is imposed, violations of LFU 
universality in K and ir semileptonic decays turn out to be unobservables, given the additional 
suppression factor \V u b(^l q )*/V uq \, for q = s,d, compared to Eq. (19). 

The universal 30% excess to r semi-leptonic charged-current decays is likely to explain, at 
least in part, the tension between exclusive and inclusive determinations of \V c b\ and \ V u b\. The 
argument goes as follows: the B —> X C}U tv decays followed by r — > Xivv represent a background 
for the inclusive B —> X c ^v analyses. At present, this background is subtracted via montecarlo 
simulations that assume a SM-like B(B — > X c ^ u tv) [32]. This procedure therefore underestimates 
the background events and leads to an enhanced B —> X c ^ u lv signal. On the other hand, the 
problem is not present in the exclusive decays of the type B —> where the kinematical 

closure of the events prevents the contamination from r decays. A precise estimate of the effect 
would require a re-analysis of B —> X c ^ u iu data and is beyond the scope of the present paper. 
However, we note that this effect necessarily goes in the right direction (enhanced signal in 
inclusive modes), and that is likely to be larger in b —> u compared to b — > c, given the different 
kinematical cuts. We are then led to the conclusion that the most reliable estimates of \V c b\ and 
| V u b\ are those obtained by means of exclusive decays and, more specifically, exclusive decays 
into electron final states. 

As a result of this discussion, we urge the experimental collaborations to reanalyze all semi¬ 
leptonic charged-current B decays without imposing LFU, both as far as signal and as far as 
background are concerned. 


3.2 Bounds from A F = 2 


Also in the case of A F = 2 transitions the new interaction amounts to an overall flavor non- 
universal rescaling of the SM amplitudes. It is therefore convenient to define the ratios 


pAF=2 

Ad 

n q 


A(B q —> B q ) SM+NP 
A(Bq —> Bq) SM 


— 1 + 


9q 




bq’ 


91 (VtlVtq) 


x 


SM 


( 21 ) 


where 5 

4°° p = a 4 5 ° 2 (Xt) « 6.5 x lO " 3 . (22) 

The consistency with experimental results on down-type A F = 2 amplitudes, where no signifi¬ 
cant deviations from the SM are observed (up to the 10%-30% level depending on the specific 
amplitude) implies |A®-| ^ 10 —1 1V 3 y|, for Rq = 0.15 and gi/g q = 0(1). As anticipated, this 
justifies the expansion on the r.h.s. of Eq. (17). 

If the corrections of A|- from the leading term are generated by U(2) q breaking spurions, as 
proposed in Eq. ( 6 ), the R^ =2 terms should respect the U( 2 ) 3 prediction [29] 

= 2 = = 2 (23) 

5 For the SM amplitude and the definition for the loop function So(xt) = Sq(m %/« 2.4 see e.g. Ref [35]. 
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that, in turn, implies no corrections to the clean ratio Furthermore, if the ei 

parameter in the expansion ( 6 ) is real, we expect no corrections to the CP-violating phases of 
B s and B<i mixing. 

To discuss the NP impact in b —> s£ + l~~ decays we need to establish a precise upper bound on 
|A{J S |. According to the U(2) 3 fit of meson-antimeson mixing of Ref. [34], we set Rg F=2 € [0.8,1.2] 
at 95%CL, that implies 


l A LI < l A LI max — 0.093 \Vt s \ 


9 A 1/2 /0T5V 72 

Sq \ RO ) 


(24) 


Even in the limit of negligible A?- for i / 3 or j / 3, a potentially sizable contribution to 
AC = 2 is generated by CKM mixing, starting from the leading term in Eq. ( 6 ). This can be 
written as 


AC 


0 v ) 

AC=2 


Rue 


* \2 


1 {VvbV£) 

r* 
cb 


A 2 

- LI -UC 


V ubV* 12 


^R 0 9 -l\V ub V: b \ 2 
LV2 si 


{uLlfiCL) 2 +h.C. , 
1 - 1/2 


6.9 x lO ' 3 TeV x 


1/2 


0.15' 


1/2 


Ro J 


(25) 

(26) 


Remarkably, for gi/g q = 0(1) this result is compatible with the existing bounds from CP 
violation in D-D mixing that require R uc > 3 x 10 3 TeV [36]. For Rq = 0.15 this fixes ||^| < 5.4. 


3.3 Bounds from LFU and LFV in r decays, and neutrino physics 

LFU in r decays has been tested at the permil level. Assuming A- is negligible if i = e or j = e, 
and imposing |A?| <C I 3 * Vy-,-], such tests can be used to set stringent limits on |A^J and |A^|. 
Moreover, a strong limit on the product |A^||A^| follows from the upper bound on B{r —> 3/i). 
The relevant modified effective Lagrangians are 


AC m ~ 

^-^LFU — 

s 2 

2 rriy 

> 

h> 

< 

11 

1 ( 

RV 


- hKn,\ 2 ) + ^Af M (r L 7 At /i L )(i/ T 7 M z/ T ) + h.c. 


j2~(T L 'y ll g L )(g L 'y ll g L ) 


R-Tfl — 


33 R \l 

IV 2 % q r r 


1 - 1/2 


(27) 


As far as LFU tests are concerned, the observable we consider is [37] 


B[t -t- guu)f(x 2 e ) 
B(t -> evu)f(xl) 


I + Rqj f x i 


- -I \ l I 

/Z/X 2 ' 


+ 


Ro Si 


= 1.0040 ± 0.0032 


(28) 


where xi = f (x 2 ^) //(xg) = 0.9726 is a phase space factor, and we summed over neutrinos of 

arbitrary flavor. The numerical result on the r.h.s. of Eq. (28) is obtained using PDG data [5]. 
Expanding to first order in Rq and assuming | A!; ;/ 1 2 <C |A^J we obtain 



= ( 0.013 ± 0 . 011 ) x (] R 


(29) 






















This constraint is significantly stronger than the bound from b —> cg(e)u universality (Eq. (20)), 
unless g q /g e > 1. 

In principle, an independent bound on |A^| can be obtained by neutrino trident production, 
namely muon pair production from scattering on a heavy nuclei. The inclusive cross section 
( t(v^N —> v fl Ng + g~) is proportional to the combination of effective couplings (Cy + C\) where 
Cy,A = Cy 1 )! + ACv,a [38], and the SM reference values are Cy M = 1/2 + 2 sin 2 6\v and = 
1/2 [38]. The corrections to these couplings in our model are 


A C v = AC A 


m W9e lyl 
2 rrfyg 2 ' w 


2 


(30) 


Combining the reported cross section measurements from the CHARM-II collaboration (a/a SM = 
1.58 ± 0.57 [39]) and the CCFR collaboration (cr/cr SM = 0.82 ± 0.28 [40]), we find 


|Aj„| < 1.5 


9t \Ro J 


(31) 


that is well below the LFU bound in Eq. (29). 

As far as LFV is concerned, the B(t -A- 3 g) < 2.1 x 1CU 8 bound [5] implies A Tfl > 11 TeV, 
that can be translated into 


|A^A^| < 0.005 


9q 

9t 



(32) 


If|AjU assumes the maximal allowed by Eq. (29), we are left with the bound |A^| ;$ 0.15. The 
latter is compatible with the hypothesis that A (A admits an expansion similar to that of AT in 
Eq. (6), or that the leading breaking of the U{ 2)^ flavor symmetry is determined by spurions 
transforming as U(2)# doublets. 


3.4 New-physics effects in b —> s£ + £ 

The effective Lagrangian encoding NP effects in b — y s£ + £~ (£ = e, ji. r) is 


AC 


(V) 

b^-s£+i~ 


2 G f 




(33) 


Using AC^^ f+f! _ to determine modified matching conditions for the Wilson coefficients of the 
most general b -A s£ + £~ effective Hamiltonian, 


-T/b-tS 


4 G F y y* e 

AM 


(Cfof + Cf'O u ) + h.c 


<£/ 


167r 2 


(34) 


leads to 


ACg = -AC[ 0 = - 


ttR 0 A? 


bs 




VtlVts 


■1M e ) _ 

y 9 


^(e) _ 
10 


' v ^(ee)^C[ 0 , 


(35) 
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where 


°i = (sL'fb L )£'y„£ , Of o = {s L Y b L)£'yui b £ ■ (36) 

Present b —> sg + anomalies and R^ e seems to indicate a LF non-universal modification 
in the Wilson coefficients Cg compared to the SM (see e.g. Ref. [13,14]). However, a good 
fit to present data is also obtained assuming ACg = -AC^ / 0 and AC| 10 = 0, that is 
compatible with the modification expected in our NP framework for |A^J | Af e |. The best fit 
thus obtained implies ACg = — AC^q = —0.53 ± 0.18 [41], 

In order to reproduce this result within our model we must impose 

= ( 3 - 4 ± 1-1) x 10“ 4 x ■ (37) 


This result is in some tension with the bounds on | Xf lfl | and | A^ s | dictated by LFU in r decays and 
A nriB q mixing, respectively. To express this tension more clearly, it is convenient to normalize 
Eq. (37) to the maximal value of |A][J allowed by A rriB q mixing. This leads to 


Ks 

( Ro N 

I V2 

9l 

Afamax 

\0-15y 


9q 



(0.09 ±0.03) , 


(38) 


that should be compared with the constraint on |A^J from Eq. (29). Given the different scaling 
of Eq. (38) and Eq. (29) in terms of gz/g q , the tension decreases for \gc/g q \ < 1. 

As far as b st + t~ decays are concerned, for Rq = 0.15 and g q = gz, we find 

AC^ = —AC [ 0 « -5.6 x , vs. (C 9 ") sm ~ -(C[ 0 ) SM ~ 4.2 . (39) 

Thus if A q bs is close to |A/ s | max , as favored by b —> sg + g~ anomalies and R^ e , depending on 
arg(A^ s ) we have two very different non-standard predictions for b —> st + t~ decays. In the 
case of maximal constructive interference of NP and SM amplitudes, b — > st + t~ rates could 
be enhanced up to a factor ~ 5 over the SM; in the case of maximal destructive interference, 
b —> st + t~ rates could be strongly suppressed (even less than 1/10) compared to the SM 
expectation. This possible enhancement or suppression would hold also for the b —> sv T v T rates, 
but it would appear “diluted” by a factor of « 3 in the measurable b —> svu rates summed over 
all neutrino species. 

In principle, the effective four-fermion Lagrangian in Eq. (4) could allow also FCNC-LFV 
transitions of the type b —>• s£±£j, with the largest amplitude expected for b —> The 

latter can be estimated by means of Eq. (35), with the replacement A^ (/J —> X e T ^. Given the 
constraint on (A/^l in Eq. (32), we find that FCNC-LFV helicity-conserving transitions (B —> 
Kt±h t , B —> it *T =t ^ =F , ...) can have rates which are at most 10% of those of the corresponding 
di-muon modes in the SM. Similarly, we find B(B S —> r ± /T F ) < 10 -8 . These bounds makes the 
experimental search of these FCNC-LFV transitions very challenging, at least in the short 
term. We also note that such bounds are saturated only if B(t —> 3/x) is just below its current 
experimental bound. 
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Obs. Oi 

Exp. bound (fit A a i) 

Def. Oi(x a ) 

Ro(D*) 
Ro(D ) 

AJCe 

0.14 ±0.04 

0.19 ±0.09 

0.00 ±0.01 


A i?|f= 2 

0.0 ±0.1 


A CS 

-0.53 ±0.18 

- (vr/a em ) X q bs /\V* b V ts \ 

ARt—> fj</e 

0.0040 ± 0.0032 

^{x^-\\xU 2 ) 

K 2 » 

(0.0 ± 4.1) x 10 -9 [GeV 2 ] 

(G F /V2)e 2 A^A^ 

K ? 

(0.0 ± 5.6) x 10“ 14 [GeV 2 ] 

(G F /V2)e 2 q \V ub V;f 


Table 1: Observables entering in the fit with their experimental bound (assuming the uncer¬ 
tainties follow the Gaussian distribution) and the expression in terms of the parameters of our 
model. 


3.5 Combined fit and discussion 


The low-energy observables discussed above depend on the three flavor-non-universal couplings 
X q s , A^ /t , A(- m , and the two flavor-independent combinations 


9e,g mw 
gm v 


122 GeV 


9t,q 


my 


(40) 


which we assume to be bounded by |e^ i9 | < 2. We have performed a combined fit of these 
parameters using the experimental constraints reported in Table 1. For simplicity, we have 
assumed Gaussian errors for all the observables. The preferred region of the model parameters 
(x a ) has been determined minimizing the x 2 distribution 


X 2 (x a ) 


E 


(Oi(x Q ) - m) 2 

■2 


(41) 


The best-fit point is found for 

e e « 0.37 , e q « 0.38 , X q bs « 2.3 x 10" 3 , X f j jp « 2.0 x 10~ 2 , X l Tfl « 4.8 X 10 -2 . (42) 


The x 2 improvement of the best-fit point with respect to the SM limit is x 2 (-' c sm) — X 2 ( x bf) = 
18.6 for 5 d.o.f., which corresponds to a p-value for the SM hypothesis of 0.002. In Fig. 1 we 
show the 68%CL and 95%CL regions in the (e q ,eg), (X q s , Aj) ;i ). (Xj lfJt , A(- ; J, and ( ACg,AR Fs ) 
planes, after having marginalised over the other parameters. 

The best-fit point implies a small non-standard contribution to Cg. This is because of 
the bounds on |A^ /t | and |A®J dictated by LFU in r decays and A ms q mixing (see sect. 3.4). 
However, in the 95%CL (68%CL) preferred region of the model parameters the effective coupling 
i±i can exceed 0.10 (0.05). In this case ACg can be within la or 2 a of its central value (see 
right panels in Fig. l). l> 

6 A “perfect fit” of ACg can be obtained extending the minimal version of the model, at the cost of introducing 
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Figure 1: Results of the low-energy fit in Table 1: 68%CL (green) and 95%CL (yellow) allowed regions 
in the (e 9 , q) plane (upper-left plot), (A)( s , ) plane (upper-right plot), (A* , ) plane (lower-left plot), 

and in the (ACg , ARb b ) plane (lower-right plot), after having marginalised over the variables not shown. 
The black dots represent the best-fit points for these 2d likelihoods. In the upper-right plot, the solid, 
dashed, and dotted red lines represent the iso-lines respectively for the best-fit, 1- and 2-cr ranges for 
ACg , with fixed Rq = 0.15. In the lower-right plot, the dashed and solid blue lines represent the 68%CL 
and 98%CL regions for ACg and A R Bs as favored by b — ► s/r + ^ and Aitib s data. 
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In summary, we find that the effective Lagrangian in Eq. (4) provides a significantly improved 
fit to low-energy data. It is worth stressing that, even from the EFT point of view, this model 
is highly constrained given the underlying set of dynamical hypotheses. As a result, the model 
leads to a series of predictions often different (more precise) than those obtained using more 
general EFT approaches (see e.g. [18,19,28]). The main predictions, which can be used to test 
the model in a more stringent way in view of future data, can be summarized as follows: 

Charged currents. The b —> c(u)tv charged currents should exhibit a universal enhancement 
(independent of the hadronic final state). This implies, in particular, Rbtu = R^ = Rgjt- 
LFU violations between b —> c{u)^v and b —> c(u)eu can be as large as 0(1%). The 
inclusive V%| and \V u b\ determinations are enhanced over the exclusive ones because of 
the r contamination in the corresponding samples. 

FCNC. The modification of the b —> si + t~ operators are purely left-handed. This implies, 
in particular, ACg = —AOfg, hence a suppression (in the 10%-20% range) of B{B St d —> 
fi ~) rates compared to their SM expectations. The central value of the anomaly in R^ e 
is likely to decrease (to ~ 10% or less). The NP contribution to the b —> st + t~ amplitude 
is likely to be close, in magnitude, to the SM one. This implies a rate enhancement of at 
most a factor of ~ 5 compared to the SM (constructive interference) or a strong suppression 
(destructive interference). The magnitude of the FCNC-LFV transitions b —> is at 

most 10% (in the rates) compared to the b —> s£ + l~ ones. 

Meson-antimeson mixing. A 0(10%) deviation from the SM is expected B s mixing, if the 
anomaly in R^ c persists. According to the most plausible breaking structure of the U(2) q 
symmetry, this deviation should be present also in Bd mixing and should preserve the 
relation AMb s /AM s d = (AMb s /AMB d )sM- The D-D mixing amplitude should acquire 
a CP-violating phase, whose magnitude could be just below the current experimental 
bounds. 

t decays. The r —> 3/r and ^ —> 3e processes are generated at the tree level (contrary to 
LFV dipole transitions l 7i —> ij'y) and could be close to the present experimental bounds, 
although no precise correlations with other observables can be derived at present. If the 
anomaly in R^ e persists, violations of LFU in r -> n&v vs. r —> eVv are expected to be 
just below the current experimental bounds. 


4 Constraints on the dynamical model 

4.1 Bounds from LEP-I 

Since the couplings of the heavy vector with SM fermions in this model are strongly non- 
universal, we cannot apply the LEP-I constraints as encoded in the bound on the S-parameter. 
Instead, we consider the non-universal fit of LEP-I data performed in the context of dimension-6 

more free parameters. In particular, a natural extension is obtained with the inclusion of a SU(2)l singlet, coupled 
to the current obtained by J£ in Eq. (5) with the replacement T a —> l,g q t —> g' q e- 
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Figure 2: Results of the combined flavor and electroweak fit: 68%CL (green) and 95%CL (yellow) 
allowed regions in the (e g , e#) plane after having marginalized over the other parameters. The black dots 
represent the best-fit points. 


operators in Ref. [43]. To do this, we translate the effective Lagrangian of Eq. (13) to the Higgs 
basis used for the fit. 

Since the constraints on the Z couplings to third generation (left-handed) quarks and leptons 
are of the same order as the bounds on the couplings to lighter fermions, while in our model 
the third generation is the one with biggest couplings, the strongest constraints on the model 
will arise from the bounds on Z couplings to third generation fermions. This motivates us to 
simplify the analysis of the EFT fit by neglecting XfJ for i. j ^ 3. In this limit the fit only 
depends on these two combinations of parameters: 


Q£H = 


g 2 rriy 


(4.3 ±8.7) x 1(T 4 , 


£h 


9qdH rn ^Y 

g' 2 my 


(-0.8 ± 1.4) x 1(T 3 , (43) 


and the correlation is negligible. We introduced the adimensional parameters ex = g.xmw / gmy , 
with X = £,q, H. With this notation the constraint in Eq. (18) from charged current R-decays 
can be written as Ro = e^e q = 0.14 ± 0.04. In Fig. 2 we combine these experimental constraints 
with the ones from flavor physics and show the 68%CL and 95%CL allowed regions in the (e 9 , en) 
plane. From this we conclude that \en\ < 5 x ICE 3 . This result allows us to conclude that, in 
absence of new degrees of freedom in the model, the massive vectors decay dominantly to SM 
fermions by means of the last term in Eq. (12). 


4.2 High-energy searches 


We parametrize massive vector boson couplings to SM fermions (in their mass-eigenstate basis) 
as follows 


ACyj = V“J“ = cYj fin*flVp . 


(44) 
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With this definition, the two body V —> fifj decay width is 


nv ->■ fifj ) 




f ™ fi ; m fi 

\my ’ my 


(45) 


where 


,) = (i - ^ V^W^TW- V -w . 


Nq is the dimension of the color representation of the fermions, and we have assumed 
m fi +m fj . 

Due to the approximate U(2) q x U(2)i symmetry, the total decay width of the vector 
is dominated by decays to third-generation fermions. In the limit rriy 2> 4ru|, 


(46) 

my > 
bosons 


m v ± 


r yO 
myo 


1 

487T 


(5l + 3ff Q 2 ). 


(47) 


The neutral vector boson predominantly decays to t + t~, v t v t , bb, and it final states. The 
relative impact of the leptonic and hadronic decay modes is driven by the ratio gi/g q (and the 
phase-space corrections to the large top-quark mass). The decay to a muon pair is parametrically 
suppressed by the smallness of . In particular, the following relation holds 

B(V° -> g + g~) = |A£/ B(V° -> r+r") . (48) 

The dominant charged vector decay modes are and t + v t . In the following we assume my+ > 
mt, such that the V + state cannot be produced on-shell from top decays. We checked by explicit 
computation that when this criteria is satisfied the corrections to t —> bT + u r decay are well below 
present experimental sensitivity. The decays of both charged and neutral states to SM gauge 
bosons are strongly suppressed due to the strong limits on the eu parameter from electroweak 
precision data (see Sect. 4.1). 

The single vector bosons production at the LHC is dominated by Drell-Yan type processes, 
i.e. pp —> V + X, where X stands for additional hadronic activity. While resonance searches 
in general impose severe limits on sequential (SM-like) W' and Z' bosons, we find significantly 
milder limits within our model. This is because of the specific flavor structure which suppresses 
both the production cross section and the decays into muon or electron final states. 

In order to derive the present collider limits on the model, we have confronted the predictions 
of the model to a number of ATLAS and CMS searches for heavy W' and Z' resonances [44-50] . 
To this purpose, we have implemented the model in Universal FeynRules Output (UFO) [53] 
using the FeynRules [52] package version 2.3.1. We have used the MG5_aMC_v2.2.3 [54] package 
to simulate the tree-level pp —> and pp —> V° production at \fs = 8 TeV in the 5-flavor 

scheme. Finally, we have validated the implementation of the heavy vector triplet couplings to 
fermions by simulating decays and comparing the numerical results with the analytic expressions 
in Eq. (45). 

In Fig. 3 (left) we show the predicted cross sections (in pb) for pp —> (pp —> V + plus 

pp —> V~) and pp —> U°, as obtained in the limit A?- = (i.e. exact flavor symmetry, but 

for the breaking terms induced by the SM Yukawa couplings). In this limit the production cross 
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sections are completely determined by g q and my. As can be seen, the neutral cross-section 
is about 100 times larger than the charged one. This is because the V° state is produced by 
bottom-bottom fusion, that is allowed in the limit of exact flavor symmetry, while the leading 
production channel is bottom-charm fusion, that is suppressed by \Vcb\ at the amplitude 
level. 

The search for W' — > tb with the ATLAS detector at 8 TeV and 20.3 fb -1 of data excludes 
a left-handed W' boson with a mass of 500 (1000) GeV if a(pp —> W') x B(W' — > tb) > 
3.3 (0.19) pb [44], In addition, the CMS search for W' —> tv performed at 8 TeV with an 
integrated luminosity of 19.7 fb -1 excludes a W' boson with mass 300 (500) GeV if a(pp —> 
W') x B(W' — > tv) > 4 (0.1) pb [45]. Comparing these limits with the predicted cross sections 
shown in Fig. 3 (left), we conclude that these searches have little impact on our model. We also 
checked that the ATLAS search for W' —t pv [46] has no relevance due to the limit on the | A^ (// 
coupling. 

The resonance searches for neutral vector bosons are more relevant due to the larger expected 
cross section. The ATLAS search for Z' resonances decaying to t + t~ using 19.5 — 20.3 fb -1 of 
8 TeV data [47] sets an important constraint on the parameter space of the model. The exclusion 
limits (under the assumption of a narrow resonance) are shown in Fig. 4 in cyan solid (dashed) 
line assuming B{Z' —t t + t~) = 0.01 (0.10). The region above these lines is excluded at 95%CL. 
The exclusion limits start from my = 500 GeV because Ref. [47] reports the limits on a x B 
only above this mass. To overcome this problem, we urge the experimental collaborations to 
extend the search for Z' resonances even in the low mass region. To extract the present limits 
for my < 500 GeV, we use the CMS search for the neutral MSSM Higgs boson decaying to a 
pair of tau leptons [51] at 8 TeV and 19.7 fb -1 luminosity. The collaboration reports a model 
independent limit on the 66-induced production cross section times B(H —» t + t~) (assuming 
a narrow resonance) in the region 100 GeV < m h < 1 TeV. We have performed a parton- 
level MadGraph simulation to compare the kinematics of the t + t~ final state produced by a 
scalar and a vector resonance (of mass 200 GeV). Having found small differences, we have re¬ 
interpreted the CMS bound into the er X B limit for our model reported in Fig. 4. In particular, 
the region above pink solid (dashed) line is excluded assuming B(Z' —y t + t~) = 0.01 (0.10). 
While the searches for dimuon resonances are usually more sensitive than the r + r~ ones (see for 
instance [48]), we find them less relevant for our model due to smallness of X^ . On the other 
hand, dijet [49] and tt [50] resonance searches set limits on the cross section times branching 
ratio for my ~ 1 TeV of the order of 1 pb. 

The impact of the direct searches in the t + t~ channel on the parameter region preferred 
by flavor data is illustrated in Fig. 4. The wide light green region is obtained imposing Rq = 
0.14 ± 0.04 (68% CL region) and g q ,gg < The narrower dark green (yellow) band is the 

region for which Rq is within 68% CL (95% CL) and g q = gg. In the minimal model, the 
predicted V° —t t + t~ branching ratio for g q = gg is B(V° —> t + t~) ~ 1/8. Comparing with 
the exclusion curve obtained for B(Z' —t t + t~) = 0.1 we deduce that the minimal model is 
ruled-out for g q = gp. The situation improves for gg <C g q (a configuration also preferred by 
flavor data, albeit with the lower bound gg > g q /5A from D-D mixing, see sect. 3.2), given 
B(V° — > t + t~) ~ (1/8) x (ge/g q ) 2 - However, it is not possible to completely evade the bounds 
for perturbative values of the couplings (g q ,ge < v/br) in the region preferred by flavor data. 

There are different ways to evade the LHC limits on pp —> V° —t t + t~ going beyond 
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Figure 3: Predicted cross sections for single and pair production of charged and neutral vector 
bosons for 8 TeV pp collisions as a function of g q and my in the limit of approximate flavor 
symmetry. 



m v o (GeV) 


Figure 4: 


Preferred region from flavor data and exclusion limits from LHC. See text for details. 
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the minimal model. The simplest possibility is to add new V° decay channels, say to a dark 
sector. This would result into lower values of B(V° — > As can be seen in Fig. 4 if, for 

9q = 9ii B(V° —> t + t ~) decreases to 0.01, then there are regions of the parameter space that are 
allowed, both at low and at high masses. Another option is to consider a heavy V° in the limit 
of a strongly coupled theory (Ty ~ My). In this case the resonance becomes broad and the 
limits obtained assuming a narrow state no longer holds. A third possibility would be to add an 
additional neutral heavy vector, e.g. a SU(2)l singlet, close in mass to the neutral component 
of the triplet, with couplings tuned to interfere destructively with U° in the pp t + t~ + X 
cross section. 

We finally comment about the pp -» VV (pair production) process at the LHC. This proceeds 
via: i) f-channel quark-exchange diagrams, controlled by Eq. (44), and ii) s-channel diagrams 
with off-shell SM electroweak gauge bosons. In the limit of no mixing with EW gauge bosons, 
9h ~ 0, and neglecting contributions from additional non-minimal operators [42], the relevant 
interactions of the heavy vectors with the SM electroweak gauge bosons described by Eq. (12) 
are 


C D -ig( 8 e A» + c e ZjV-{d*V+ v -PV+^-igW+tV^^V- -d v V~) 

MigW+V -- d u V 0fl ) + h.c. , (49) 

where g is the SU{2)l coupling constant and sg ( eg ) is sine (cosine) of the Weinberg angle. As 
for the single V production, we have used MadGraph to simulate pp -» V + V~ and pp — > V^V 0 
(that is, pp —> V + V° plus pp —> V~V°). In Fig. 3 (right) we show the predicted cross sections 
(in pb ) for 8 TeV proton-proton collisions, as obtained in the limit X')- = Sisd^j. As illustrated 
with the vertical iso-lines, the cross sections are dominated by s-channel diagrams for small g q 
couplings. On the other hand, for large couplings there is a substantial contribution from the 
bb induced f-channel diagram to V + V~ production. Based on the predicted cross sections, we 
conclude that the single production is more relevant compared to the pair production for the 
interesting region of the parameter space. 

5 Conclusions 

Lepton Flavor Universality is not a fundamental symmetry: within the Standard Model it is an 
approximate accidental symmetry broken only by the Yukawa interactions. This specific sym¬ 
metry and symmetry-breaking pattern results in tiny deviations from LFU in helicity-conserving 
amplitudes, within the SM, and it implies that LFU tests are clean probes of physics beyond 
the SM. 

Motivated by a series of recent experimental results in B physics pointing to possible viola¬ 
tions of LFU, both in charged and in neutral currents, in this paper we have consider a simplified 
dynamical model able to describe these effects in a unified way. In particular, we have shown that 
a SU(2)l triplet of massive vector bosons, coupled predominantly to third generation fermions 
(both quarks and leptons), can significantly improve the description of present data. 

The proposed model has a series of virtues compared to previous attempts to describe such 
effects in terms of New Physics: i) it connects the breaking of LFU between charged and neutral 
currents, and between semi-leptonic and purely leptonic processes; ii) it is based on a simple 
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flavor symmetry, whose breaking terms are related to the structure of the SM Yukawa couplings, 
both in the quark and in the lepton sector; iii) it connects low-energy deviations from the SM 
to direct searches for NP at high px- The constrained structure of the model makes it highly 
non trivial to satisfy all existing bounds and, at the same time, accommodate deviations from 
the SM as large as indicated by the central values in Eqs. (l)-(3). We find that this happens 
quite naturally in the case of charged currents, both at low and at high energies. The situation 
is more problematic in the case of neutral currents. On the one hand, the maximal deviations 
from unity in R^ e barely exceed 10%. On the other hand, the minimal version of the model 
is ruled out by the direct searches for resonances decaying into t + t~ at ATLAS and CMS. As 
discussed, both these issues can be improved with less minimal versions of the model, at the 
cost of introducing more free parameters. 

One of the most remarkable aspects of the minimal version of the model is the well-defined 
pattern of deviations in low-energy processes listed at the end of Sect. 3.5. This pattern is largely 
insensitive to possible extensions of the model necessary to overcome the constraints from direct 
searches. It mainly reflects the symmetry structure of the model and could be used, in the near 
future, to verify or falsify this framework with more precise data. 

Besides the specific predictions of the proposed model, we stress the importance of future 
experimental tests of LFU at low-energies, and dedicated searches for flavor-non-universal phe¬ 
nomena at high energies. On the low-energy side, we urge the experimental collaborations to 
re-analyze charged-current B decays without assuming lepton flavor universality. On the high- 
energy side, we encourage the search for deviations from the SM in t + t~ and tt invariant-mass 
distributions, relaxing the hypothesis of narrow resonances and covering also the region of low 
invariant masses. 
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